Successive phase transitions at finite temperatures of the supersolid in the 
three-dimensional extended Bose-Hubbard model 
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We study the finite temperature properties of the extended Bose-Hubbard model on a cubic lattice. 
This model exhibits the so-called supersolid state. To start with, we investigate ordering processes by 
quantum Monte Carlo simulations, and find successive superfluid and solid phase transitions. There, 
we find that the two order parameters compete with each other. We obtain the finite temperature 
phase diagram, which contains the superfluid, the solid, the supersolid and the disordered phase. 
We develop a mean-field theory to analyze the ordering processes and compare the result with that 
obtained by simulations, and discuss the mechanism of the competition of these two orders. We also 
study how the supersolid region shrinks as the on-site repulsion becomes strong. 



INTRODUCTION 

The supersolid state is an interesting state of matter 
which has both solid and superfluid properties. The solid 
state is characterized by the breaking of translational 
symmetry, and the superfluid state is characterized by 
the breaking of U(l) symmetry of the phase of macro- 
scopic wave function. Thus, the simultaneous breaking 
of these two symmetries indicates that there is a flow 
component in solid. The possibility of the supersolid was 
first discussed by Penrose and Onsager [1|. Since then, 
various studies on supersolid have been conducted from 
both experimental and theoretical points of view. 

As to observation of supersolid, Leggett suggested that 
non-classical rotational inertia (NCRI) would be avail- 
able to detect supersolid in rotating solid 4 He. Recently, 
Kim and Chan [2J], reported that they found NCRI in 
solid 4 He. Although it has been pointed out that the ob- 
served NCRI may not due to the supersolid but due to 
the grain boundaries between polycrystals [3(, the topic, 
however, still attracts researchers' interest. 

The possibility of supersolid on lattice models has been 
discussed actively. Andreev and Lifshitz [J], suggested 
that the derealization of the vacancies in crystal causes a 
mass flow. Matsuda and Tsuneto [5(, studied the ground 
state of the hardcore Bose-Hubbard model using mean- 
field theory, and they showed that the supersolid is possi- 
ble when the interaction of particles has frustration. This 
is confirmed by numerical simulations in two-dimensional 
[a, ll| , and three dimensional cases [8| . 

It has been pointed out that the supersolid state can 
be realized even in non-frustrated lattice if double oc- 
cupancy of the particles is allowed [9( . In the case of a 
square lattice, the supersolid is found in the ground state 
when the number density of particle is more than that 
of the half-filled case (p > 1/2) [§]. Although ground 
state properties of the supersolid have been studied ex- 



tensively, there is no direct study on the finite temper- 
ature properties of supersolid problem in the extended 
Hubbard model. 

In this paper, we study finite temperature properties of 
the supersolid state in the extended Bose-Hubbard model 
on a cubic lattice. In our model, the hopping term causes 
the superfluid order and the nearest neighbor repulsive 
interaction tends to form the solid order. These orders 
compete with each other, and in the hardcore limit, these 
orders cannot be realized simultaneously. Thus, there 
is no supersolid state. However, as has been pointed 
out [9(, these orders can be realized together in the soft- 
core case even in not-frustrated lattices such as a square 
lattice. We investigate ordering processes of superfluid 
and solid at finite temperatures in a cubic lattice, us- 
ing the stochastic series expansion (SSE) method. As 
a result, we have a phase diagram of superfluid, normal 
solid, disordered state, and supersolid. We also study the 
finite temperature dependence of the orders by making 
use of mean-field (MF) analysis and compare the result 
with that obtained by SSE. They show qualitatively good 
agreement. Moreover, the competition of the solid and 
superfluid order is discussed using Ginzburg-Landau free 
energy. Finally, we study the effect of on-site repulsion 
U on the coexistence of the two orders. We find how 
the supersolid region at finite temperature shrinks as U 
becomes large. 



MODEL 

We analyze the extended Bose-Hubbard Hamiltonian 
on a cubic lattice 

H = —t 2J ( a l a j + (Ha}) + V 22 niU i 



~zU^2ni(ni- l)-/jt^2', 



(1) 



where aj and a^ are the creation and annihilation oper- 



ators of a boson ([a 



,tl 



Sij), and rii 



a]a,i. 



The 



parameter t denotes the hopping matrix element, U and 
V are the on-site and nearest neighbor repulsion, respec- 
tively, and [i is the chemical potential. The notation (ij) 
means the sum over the nearest neighbor pairs. The sys- 
tem size is N = L 3 , where L is the length of the system. 
The order parameter of the solid state is 
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where Q = (n, 7r, 7r) is the wave vector that represents the 
staggered order. The order parameter of the supcrfluid 

state is 



Ps = 



-^^( a l a k+ a i a j)> 



N 2 



(3) 



jk 



which represents the off diagonal long range order 
(ODLRO). 

METHODS 

We use the following two different methods to analyze 
properties of the system. 

(1) Stochastic Series Expansion 

We preform numerical simulation of the stochastic se- 
ries exp ansion (SSE), which was invented by Sandvik 
[id lll| . This method is one of quantum Monte Carlo 
(QMC) simulations, and has been successfully applied 
for various quantum systems. In order to avoid the clus- 
terization due to diagonal frustration, we adopt the gen- 
eralized directed loop algorithm [12J. We use a package 
of the ALPS [l3j. We adopt a simple cubic lattice of 
N = L 3 sites with periodic boundary conditions along 
all the lattice axes. In the simulation, we estimate the 
superfluidity p s , Eq. ([3]), using the winding number W 
of world lines defined by, 



„SSE 
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(4) 



which represents well p s 
temperature. 



|l4 , lla ] . Here (3 is the inverse 



(2) Mean-Field Approximation 

We also analyze the ordering processes by the mean- 
field (MF) approximation [16[. In order to study the 



solid state, we use a sublattice structure which is char- 
acterized by a staggered order of the density. Here we 
adopt mean-fields for the solid order and superfluid or- 
der at sublattices A and B. The Hamiltonian for this MF 
is given by 



(mf 


= Ha + Hb + C, 


(5) 


Ha 
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Hb 


= —zt(a} B + a B )(/>A + zVubtua 






u , 
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C 


= 2zt(f)A<t>B — zVrriATnB, 


(8) 



where z — 6 is the number of nearest neighbor sites. 
Here, ttia and ttt-b are the mean fields correspond to the 
expectation values of the number operators for A and B 
sites, respectively Eq. ©. Similarly, 0a and 0b corre- 
spond to the expectation values of the annihilation oper- 
ators for A and B sites, respectively. 



m A = (tia), ran = («-b), 

4>A = (oa), 0B = (ob)- 



(9) 
(10) 



Ha (Hb) is a mean-field Hamiltonian at a site of the A 
(B) sublattice. C is a correction term compensating the 
double counting of the energy. The partition function 
and the free energy are given by 



Fmf = — TjhxZMF 



(11) 

(12) 



In MF, m = (t?ia — ms)/2 denotes the order parameter 



of solid and fulfill the relation, SV 



Similarly, 



(4>A + 0b)/2 represents the order parameter of superfluid 
and fulfill the relation, p s = 2<j) 2 . 



GROUND STATE PROPERTIES 

Before analyzing properties at finite temperatures, let 
us briefly summarize the ground state properties. As has 
been reported, the system has the supersolid phase in the 
ground state when U takes a finite value @, [l7| . In Fig. 
[TJ we show the ground state phase diagram in the coor- 
dinate of (t/U, n/U) obtained by MF method. Here we 
set V = U/z at which the superfluidity takes maximum 
value [18( . The solid line is the phase boundary obtained 
by MF at T/U = 0.001. This ground state phase dia- 
gram agrees well with that obtained by the Gutzwillcr 
variational method by Otterlo [17|. In Fig. [TJ we find 
four different phases, i.e. Mott-insulator (MI), normal 
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FIG. 1: The ground state phase diagram of the soft-core Hub- 
bard model for V — U/z obtained by MF. There are four 
phases: NS (normal solid), SS (supersolid) , SF (superfluid), 
and MI (Mott-insulator). Temperature dependence of order 
parameters calculated by SSE is given in Fig. [2] at the posi- 
tions denoted by circle, triangle and rectangle. 



solid (NS), superfluid (SF) and supersolid (SS). In ad- 
dition to these four phases, a disordered phase repre- 
senting the normal liquid (NL) phase appears at finite 
temperature. We also study the ground states for sev- 
eral parameter sets by SSE. For example, we find the 
solid state for the set (t/U = 0.02, fi/U = 0.7) de- 
noted by the circle in Fig. [H the superfluid state for 
(t/U = 0.08, [i/U = 0.7) by the rectangle, and the su- 
persolid state for (t/U = 0.045, /i/C/ = 0.7) by the trian- 
gle. The temperature dependences of order parameters 
on these points are given in the next section in Fig. [2] (a), 
(b) and (c). We find the dependence on the parameters 
obtained by SSE agrees with that obtained by MF. 



PHASE TRANSITIONS AT FINITE 
TEMPERATURES 

Now, we study the ordered states at finite tempera- 
tures. The phase transition between the normal liquid 
phase and the solid phase belongs to the universality class 
of the Ising model. The phase transition of superfluid 
belongs to the XY universality class. In this section, we 
study the temperature dependence of these order param- 
eters. 



Stochastic Series Expansion 

First, we show the results obtained by SSE. The simu- 
lations were performed in the grand canonical ensemble 
using a system sizes A*" = 10 3 and N — 12 3 . 

We plot the order parameters of solid, S n , and that of 
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FIG. 2: Temperature dependence of order parameters ob- 
tained by SSE (V = U/z). (a) (t/U = 0.02, (i/U = 0.7): Nor- 
mal liquid - Normal solid transition, (b) (t/U — 0.08, fi/U = 
0.7): Normal liquid - Superfluid transition. (c) (t/U = 
0.045, n/U = 0.7): Normal liquid - Normal solid transition 
and Normal solid - Supersolid transition. 



superfluid, p s , as a function of temperature for various 
values of t. In Fig. [2] (a), we show the transition from 
normal liquid to normal solid for t/U = 0.02 in which 
only Sn appears. In the same way, the transition from 
normal liquid to superfluid for t/U = 0.08 is depicted 
in Fig. [2] (b). For t/U = 0.045, the system shows suc- 
cessive transitions and the supersolid state is realized at 
low temperatures. There, we find that the solid order ap- 
pears at a higher temperature (Fig. [2] (c)). Note that the 
solid order is suppressed when the superfluid order ap- 
pears. Thus, we expected that the solid fraction and the 
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FIG. 3: The t - T phase diagram for V/U = l/z and 
fi/U = 0.7 obtained by SSE. The transition temperatures 
of the solid state Ts^ are plotted by solid circles, and those of 
the superfluid state T Ps are plotted by open circles. The lines 
connect the data points for the guide of eye. 



superfluid fraction compete with each other. It should 
be noted that p s appears at higher temperature than the 
solid order for t/U = 0.55 (not shown). 

In Fig. [3J we depict a phase diagram in the coordi- 
nate of (t/U, T/U) for the fixed values, V/U = l/z and 
[i/U = 0.7. The transition temperatures of the solid 
state Tg^ are plotted by solid circles, and those of the 
superfluid state T Pb are plotted by open circles. To de- 
termine the transition temperatures for each value of t, 
we use the method of the Binder parameter of the sys- 
tems with L=10 and 12. In Fig. [3j there are four differ- 
ent phases: NL, NS, SF, and SS. These phases meet at 
a tetra-critical point, (t c , T c ). The competition of solid 
and superfluid orders is also found in the phase diagram 
(Fig. [3]). Namely, above t c , the transition temperature 
of solid is smaller than that of smooth extension of T Ps 
from t < t c . Therefore we conclude that T Pa is suppressed 
from that of the case in which the superfluid would not 
order. Similarly, below t c , the transition temperature of 
superfluid is smaller than that of the case in which the 
solid would not order. 



Mean-Field Analysis 

Here, we calculate the temperature dependence of or- 
der parameters by making use of MF. In Fig. |4j we 
show the successive transitions of superfluid and solid 
for (t/U = 0.045, /z/C7 = 0.7). As was seen in the SSE 
simulation, here we find again the suppression of the solid 
order by the superfluid fraction. Namely, S„ has a cusp at 
the superfluid transition point. We also depict the phase 
diagram and compare that to that of SSE (Fig. [5]). They 
show a qualitatively good agreement, e.g., there is the 
tetra-critical point (£ C ,T C ), and the critical temperatures 



0.3 






Ps 






1 ? 


25 










1 


0.2 










0.8 


0.15 












0.6 


0.1 












0.4 


0.05 


"\ 




















0.2 
















.4 


3 0.1 


0.2 




0.3 









t/U 











FIG. 4: Temperature dependence of order parameters for 
V/U = l/z, t/U = 0.045, and fi/U = 0.7 obtained by 
MF. The solid line denotes the solid order parameter and 
the dashed line denotes the superfluid order parameter. 
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FIG. 5: The t - T phase diagram for V/U = l/z and 
H/U = 0.7 obtained by MF. The transition temperature of 
the solid state Ts x is denoted by the thick solid line, and that 
of the superfluid state T PS is denoted by the thick dashed line. 
We also depict the phase diagram (solid and open circles) ob- 
tained by SSE (Fig. [3} for comparison. 



, s are suppressed by appearance of the other 
as mentioned in the last session. 



T S7r and T, 
order 

Let us study the competition between the solid and su- 
perfluid orders by analyzing the Ginzburg Landau (GL) 
free energy. Since the order parameters m and <f> take 
small value in the vicinity of the tetra-critical point, the 
GL free energy is expressed as 



F = am 2 + bm 4 + c0 2 + dep* + hm 2 <j) 2 



(13) 



When a becomes zero at Tg^ ( a ~ ao(T — Tg n )) with a 
positive b, the second-order transition between the nor- 
mal solid and normal liquid phase takes place, and simi- 
larly when c becomes zero at T Ps (c ~ cq(T — T Ps )) with a 
positive d, the second-order transition between the super- 
fluid and normal liquid phase takes place. The fifth term 
represents the competition between the solid and the su- 
perfluid order. If h equals zero, the transition of the solid 
phase and that of superfluid take place independently at 
Tg„ and T Pg , respectively. 



If h is positive, the presence of the superfiuid order low- 
ers the transition temperature of solid. Below the t c , the 
solid order emerges first as temperature decreases (Tg^). 
Then, the superfiuid order appears at the modified tran- 
sition temperature T' which is smaller than the original 
one (Eq.(XD). 



K. < T P s 



a h/2bco 
1 — aoh/2bcQ 



(Ts 



T Ps )<T p 



(14) 



In the same way, the solid order lowers the transition 
temperature of superfiuid. The changes of the transition 
temperature become large when h becomes large, and 
this means the shrinkage of the supersolid region. Thus, 
h represents the competition between the solid and the 
superfiuid orders. 

The coefficients in the GL free energy can be calculated 
from the MF Hamiltonian, Eq. ([5]). For example, the 
coefficient a is expressed as 



a= - l)V - (z - 1) 2 V 2 /Z_ 



MF 



P 



xTrl d\e XHMF {n A - n B )e~ XHMF (n A - m) • (15) 

which is proportional to (T— Tg^) near the critical point. 
As a final part of this section, we discuss the effect 
of on-site repulsion on the coexistence of solid and su- 
perfiuid orders. As has been mentioned, the SS phase 
does not exist in the hardcore case. Since the hardcore 
is the limiting case of infinite U , we expect that the su- 
persolidity is suppressed when the on site repulsion U 
becomes large. We depict the t — T phase diagram for 
various values of U in Fig. [6] obtained by MF. Here, we 
use zV as the unit of energy instead of U because now 
we want to study the effect of U. In Fig. [5J we find that 
the SS region becomes narrower as U increases. Finally, 
the supersolid region disappears completely in the hard- 
core limit. Thus, we conclude U, i.e. the hardness of the 
particle, suppresses the coexistence of the two orders. 
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FIG. 6: Finite temperature phase diagrams obtained by MF 
for various values of U/V. (a) U/V = 6, (b) U/V = 10, (c) 
U/V = 100, (d) hardcore (U/V = oo). 



on-site repulsion U in particular near the tetra-critical 
point, we find that larger U enhances the competition 
between two orders and causes the shrinkage of the su- 
persolid region. It would help us to understand how the 
softness of particles contributes to the realization of the 
supersolid. 

Possibility of realization of the supersolid state on the 
optical lattice has been discussed recently. The realiza- 
tion of the Bose-Hubbard model in the optical lattice 
has been discussed [19( . For realization of the supersolid 
state, the nearest neighbor repulsive interaction V plays 
an important role. Mazzarella, et al. discussed how to 
introduce the nearest neighbor interaction [20(. We ex- 
pect that the parameters of the system can be widely 
controlled in the optical lattice and there the properties 
of the phase diagram obtained in this paper will be ob- 
served. 



DISCUSSION AND SUMMARY 



By using SSE simulation and MF analysis, we study 
successive phase transitions of the solid order and su- 
perfluidity at finite temperatures in the case the system 
has the supersolid state in the ground state. The strong 
hopping term (large t) favors the superfiuid and the su- 
perfiuid phase appears at higher temperature in the large 
t region. On the other hand, in the small t region, the 
normal solid phase appears first. In both cases we find 
that the orders of superfiuid and solid appear at different 
temperatures, i.e., successive phase transitions. We find a 
phase diagram with a tetra-critical point and studied the 
competition between the solid and superfiuid orders. By 
analyzing the dependence of the phase diagram on the 
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